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Abstract 

"Relevation"  is  the  name  given  by  Krakowski {Rev.  Franqaiee 


Autonzt.  Inform.  Recherche  Op 


given  by  Krakowski 'TI97^ 
7,  Ser.  V-2,  M-liop 


to  the  distribution  of 


failure  time  of  a  replacement  from  an  aging  stock.  The  present  authors 
(Johnson  and  Kotz  (4979t-T9BT)  {IEEE  Trane .  Reliability ,  R-28 ,  292-299; 

American  Journal  of  Mathematioal  and  Management  Soiencei  t~ If have  extended 
this  concept  to  include  (i)  hierarchal  replacement  systems  and  (ii)  dependence 
between  lifetimes  of  original  and  replacement  items.  In  this  paper,  we  present 
some  further  developments ,  including  first  steps  towards  a  synthesis  of  (i)  and 


Words  and  ftirases:  Relevation;  Dependence;  Aging  in  Stock;  Aging  in 

Service;  Farlie-Gunbel-Morgenstem  Distribution; 


Hierarchal  Systems. 


I — WSTRBUTl  ON  STAT5S3W  ~X  | 

Approrcd  *or  P'.’.M-C  rcW.se;  j 
’  Dim;?-. «. .oo.  rv-!  .--d  ! 


8  i  6 


2  5  01 


2 


1.  Introduction. 

The  concept  of  a  "relevation"  was  introduced  by  Krakowski  (1973)  to 
correspond  to  replacement  from  aging  stocks.  Suppose  that  when  an  element  with 
survival  distribution  function  (SDF)  S^(t)  fails,  it  is  replaced  by  an  element 
with  SDF  S2 (t)  taken  from  stock.  If  failure  occurs  at  time  t^,  then  the 
conditional  SDF  of  the  replacement  element  is 

S2ll(t Itl5  "  S2(t)/S2(t1)  , 

and  the  SDF  of  time  of  failure  (T)  of  the  replacement  element  is 

St(r)  -  Sj (t)  +  S2(t)  JJ  f1(t)(S2(t)}'1  dt  ,  (1) 

where  f^(t)  *  -dSj(t)/dt  is  the  probability  density  function  (PDF) 
corresponding  to  Sj (t) . 

Sj.(t)  is  called  the  relevation  of  Sj(t)  with  S2(t). 

This  concept  has  been  extended  to  situations  wherein  there  is  a  hierarchal 
system  of  replacement,  and  wherein  the  lifetimes  of  elements  are  not 
independent  (Johnson  and  Kotz  (1979,  1981)).  In  the  present  paper,  we  study- 
combinations  of  extensions  of  this  kind,  and  also  give  exanples  in  which 
differential  rates  of  aging  in  storage  and  in  service  are  allowed  for . 

2.  Dependent  Relevations  in  a  Hierarchal  System. 

Consider  a  set  of  n  elements  (Cj),  with  a  single  replacement  element  £2. 
This  might  occur,  for  example,  when  £2  called  into  service  as  a  back-up 
only  when  there  is  a  failure  among  (£j).  This  system  (comprised  of  {£j}  and 
?2)  would  quite  possibly  be  only  a  small  part  of  a  corplete  hierarchal  system; 
study  of  such  relatively  simple  systems  is  an  essential  prerequisite  to  study 
of  more  elaborate  systems  of  which  they  could  be  a  part. 


We  denote  the  failure  times  of  the  elements  of  (C^)  by  , 

and  of  by  The  joint  survival  distribution  (SDF)  of  these  (n+1) 


failure  times  is 


S12^~l;t2^  *  P^.[|1^Tlj>tlj^  n  ^■T2>t2^ 


(tj  *  (t^p. . . ,  t^)) .  The  joint  SDF  of  ^  *  minCIp, . . . ,  T^)  --the  failure 
time  of  the  first  member  of  {E^}  to  fail  --is 

P^lonTo? - ~~i 

S12^tl’t2')  =  Sl2*'tli';  l2^  I  KTrs  WA&i  7 

r"-'u:  tah  \r 


where  1*  «  (1,1,...,  1) ;  and  the  SDF  of  is 


S1  (tj)  *  S12(tjl’;  0) 


The  SDF  of  the  replacement  element  is  then 


cod 

U'icnt  ic;. 


1  *  :  ■ :  •> »  •  •  4  . 


rT  ^2  ll  |t^ 


where  f^ft)  *  -dSjftj/dt  and  S^,  j ^ Ct 1 1)  *  P(T2  >x  |T^«t) . 

In  general,  explicit  evaluation  of  (5)  leads  to  rather  conplex  formulas. 
As  a  compromise  between  realism  and  sinplicity,  we  will  (similarly  as  in 
Johnson  and  Kotz  (1981))  use  a  generalized  Farlie-Gumbel-Morgenstem 
distribution  (Johnson  and  Kotz  (1975))  as  a  special  case  for  detailed 
examination.  We  suppose 
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Si2{~l:t2 


n  n 


♦  e  Y  ?  u  -  set..)} a  -  set..,)}]  . 

jf  <  j  •  13  13 


(6) 


The  parameters  a  and  6  have  to  satisfy  certain  conditions  for  the  distribution 
to  be  a  proper  SDF.  These  conditions  can  be  suimiarized  as 


|(2r-n)aj  -  %{(2r-n)2  -  n)B  si  (r  *  0,1,.. . ,  n)  , 


or  equivalently. 


(7) 


r'|a)  -  Jjfr'2 -n)B  <  1 


r'  *  0,2,...,n  for  n  even 
r’  =  l,3,...,n  for  n  odd  , 


The  marginal  SDF  of  each  is  S(t) ;  the  marginal  SDF  of  is  S2(t).  The 
joint  SDF  of  T21,',*’Tln  *s 


Sr  CtJ  -  {  n  S(t. . ) }  [1  ♦  6  Y  l  <1  *  S(t.  •)}  (1  -  S(t..,)}]  ,  (8) 

ll  1  j*i  j  <  j  *  13 

which  is  an  exchangeable  distribution.  This  is  likely  to  be  a  reasonable 
assunption  if  there  are  no  systematic  differences  among  the  conditions  to  which 
the  n  elements  of  }  are  exposed. 

The  SDF  of  Tj  (time  to  first  failure  in  (£j))  is 

S1(t1)  -  ISCtj)}*1  [1  ♦  Hn(n- 1)6(1  -  Sftj)}2]  ,  (9) 

and  the  corresponding  density  function  (PDF)  is 


s 


*  nfCtjHSCtj)}11'1  Cl  +  *s(n -l)8(n  -  2(n  +  l)S(t1)  ♦  (n+2)  [Sftp]2}] 

-  nfCtpCSCtj)}11'1  [l  +  JiCn-DBU-sCt^Mn  -  (n^Sftj)}]  .  (10) 

The  joint  SDF  of  Tj  and  T2  is 


S12(tl»t2)  *  W{S(V}n  V  *  na{1  •  SCt^Hl  -  S2ft2)} 

+  JsnCn-l)6{l  -  S(tj)}2]  , 


(11.1) 


with  joint  PDF 


f12Ctl’V  *  nf^tp^CtpCSCtj)}"'1  [1  +  a(n  -  (n+DSCt^Hl  -  2S2(t2)} 


♦  H(n-1)B(1  -  SCtj) } (n  -  (n^Sftj)}] 


(11.2) 


The  conditional  PDF  of  T2,  given  T^,  is 


f  ft  It  )  -  f12(tl’V 

2|1° t2,tl)  '  '  fj'ctj) 


f2(t2) 


l  +  a(n-  (n+l)S(t1)}{l-2S2(t2)}+  **(n-l)B(l  -  S(tj)Hn  -  (n^jSft,)} 
1+  hCn-DBfl  -  Sftj)}  (n  -  (n^jSdj)} 


(12.1) 


and  the  conditional  SDF  is 


c  r*  e  r*  }  1  +  a{n- (n+1)SCti)Hl‘S2(t2)}+,s(n'1)B{1'S(tl)Hn'(n+^S(tl)} 

2jl  2  1  "  2  2  l*H(n-l)BU  -  SCt2)) (n  -  (n^Sfrp} - 


(12.2) 


Hence  from  (5)  the  SDF  of  time  to  failure  (T)  of  the  £2  (replacement)  element 
is 

Sj-Ct)  -  (S(r)}n  [1  +  %n(n -1)6(1  -  S(t)}2] 

♦  n  /J  f(t)(S(t)}n‘1  [1  ♦  »*(n -1)6(1  -  S(t)Hn  -  (n  +  2)S(t))] 

S2(t) 

X  — - 

s2(t) 

1  ♦  a{l  -  S2(T)}{n-(n+l)S(t1)}+  %(n-l)6(l  -  Sft^Hn  -  (n+2)S(tj)} 
1  +a(l  -  S2(t)}(n  -  (n+DSCtj))*  %(n-l)6(l  -  S(t)Hn  -(n*2)S(t)) 

(13.1) 

As  a  further  specialization,  we  will  take  S2(t)  *  S(t).  Then  making  the 
substitution  s  *  S(t)  in  the  integral  leads  to 

Sj.(t)  *  {S(r)}n  [1  ♦  Hn(n -1)6(1  -  S(t))2] 

+  nS(T)  /s(t)  s"'2!1  +  %(n-l)6(l-s)(n  -  (n  +  2)s)] 

x  1  *  c*(l  ~  S (t)  )  (n  -  (n*l)s)*»i(n-l)6(l-s){n  -  (n+2)s)  . 

1  ♦  a(l-s)(n  -  (n+l)s)  +  h(n-l)B(l-s){n  -  (n+2)s) 

(13.2) 

Note  that  Sj.(t)  is  always  the  same  function  of  S(r),  whatever  the  functional 
form  of  the  latter. 

The  integral  in  (13.2)  can  be  evaluated  in  terms  of  elementary  functions 


(see  Appendix) ,  though  the  expression  it,  in  general,  rather  complicated.  For 
the  special  case  a  ■  0  (so  that  replacement  lifetime  is  independent  of  the 
lifetimes  of  the  initial  elements). 
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Sj-Ct)  -  {S(T)}n  [1  ♦  %n(n-l)fi{l  -  S(t)}2] 


nS(T)  4)  sIl'2{1  *  Wn-l)B(l-s){n  -  (n+2)s)]4s 


-  S(T)[n-{S(T)}n-1] 


n-1 


4  l  -  \  {SCi)}0'1  [n(n  +  l)  -  2(n2  -I)S(t)  +'  n(n-l){S(r)}2J] 

n+i 


♦  B  6  . 


\  +  Bn 


(14) 


Some  numerical  values  of  A  and  B  are  shown  in  Table  1. 

n  n 


TABLE  1 

Values  of  A^  and  B^  for  S^Ct)  *  Ar  +  Bn  B  when  a  *  0. 


n 


S(T) 


\ 


n 


B_ 


n  n 


\» 


B. 


n 


n 


n 


3 


5 


This  gives  some  qualitative  appreciation  of  the  effect  of  non-zero  a. 

To  see  how  these  results  can  be  used  in  dealing  with  more  complicated 
systems,  consider  the  subsystem  set  out  schematically  in  Figure  1.  The 


I - 1  I - 1 

FIGURE  1 

component  £3  is  related  to  the  two  £2’s  ^  tlie  same  way  85  eacl1  ^2  is 
related  to  the  two  Cj’s  belonging  to  it.  comes  into  operation  at  the 
first  failure  of  the  two  £ ^ *s-  The  relevation  SDF  £p(x)  given  by  (13.1)  now 
plays  the  role  previously  played  by  the  SDF  of  each  (S(t)).  Of  course,  if 
the  resulting  formula  for  the  SDF  of  the  £3  replacement  were  to  be  written  out 
explicitly,  it  would  be  very  lengthy  indeed.  However,  the  computation  of 
nunerical  values  is  quite  straightforward.  It  can  be  carried  out  in  stages: 
first  finding  S^(x)  and  then  combining  it  with  the  SDF--Sj(t),  say --of  the 
$3  components.  (The  case  of  independent  conponents  has  been  treated  by 
Johnson  and  Kotz  (1979).) 

3.  Differential  Rates  of  Aging  in  Storage  and  in  Service. 

It  is  likely  that  an  element  will  be  under  less  stress  while  in  storage 
than  in  service.  This  can  be  allowed  for  by  regarding  a  period  t  in  storage  as 


^3 

h 


equivalent  to  a  period  g(t)  in  service  (usually  with  g(t)  <  t).  Equation  (1) 
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is  then  replaced  by 


Sj-Ct) 


Sj(T) 


S2(g(t)+  T  -  t) 

s7(g(t)) 


dt 


CIS) 


and  (2)  by 


&J.CT) 


Sj(T) 


rT  f 

1 


S2h(g(t)  +  T-t|t) 

ft)  -/l1 - - - , -  dt 

1  J  S2|l(g(t)|t) 


(16) 


One  might  reasonably  take  gft)  *  yt  --so  that  1  year  in  storage  is 
equivalent  to  y  years  in  service.  Usually  one  would  take  y  <  1,  but  this  is 
not  essential.  Even  if  it  is  assumed  that  S^t)  *  S2(t)  *  S(t),  it  is  no 
longer  true  that  &p(T)  depends  only  on  S(t)  .  As  an  example,  we  take  y  <  1 
and 


S(t)  =  1  -  ect  (OstSTS  k'1)  . 
Then  from  (15) , 

Vx) .  i  -  « .  jg  ■ 1 ^ 
*  1  -  tcry  *  -  (1  -  icry)y 


logfl  -  <Ty)  . 


(17) 
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Some  numerical  values  are  given  in  Table  2. 


TABLE  2 

Values  of  S^t)  when  S(x)  *  1  -  kt  has  specified  values 
(Y  *  service  years  equivalent  to  one  year  in  storage) . 


S(T)V 

0.2 

0.4 

0.6 

0.8 

1.0 

0.2 

0.661 

0.639 

0.611 

0.562(5) 

0.522 

0.4 

0.812 

0.804 

0.794 

0.781 

0.767 

0.6 

0.918 

0.915 

0.913 

0.910 

0.906 

0.8 

0.980 

0.979 

0.979 

0.979 

0.979 

(The  last  column  corresponds  to  no  difference  between  aging  in  storage  and  in 
service.)  The  effect  of  differential  aging  is  not  very  great,  especially  at 
the  longer  durations. 

Acknowledgement  s 

Norman  L.  Johnson's  work  was  supported  by  the  National  Science  Foundation 
under  Grant  MCS-8021704.  Samuel  Kotz's  work  was  supported  by  the  U.S.  Office 
of  Naval  Research  under  Contract  N00014-81-K-0301. 


References 

Johnson,  N.L.  and  Kotz,  S.  (1975).  Cta  some  generalized  Farlie-Gumbel- 
Morgenstem  distributions.  Com.  Statist.,  4,  415-427. 

Johnson,  N.L.  and  Kotz,  S.  (1979).  Models  of  hierarchal  replacement.  IEEE 
Trane.  Reliability ,  R-28 ,  292-299. 

Johnson,  N.L.  and  Kotz,  S.  (1981).  Dependent  nelevaticns:  Time-to-failune 

under  dependence.  American  Journal  of  Mathematical  and  Management  Sciences, 

Krakowski,  M.  (1973).  The  relevatian  transform  and  a  generalization  of  the 
gamna  distribution  function.  Rev.  Franqaiee  Automat.  Inform.  Recherche 
Op. ,  7,  Ser.  V-2,  108-120. 


11 


The  integrand  in  (13.2)  can  be  expressed  in  the  form 


a((s-b  )2  ♦  c)]'1  ?  d  ,s: 

n  n  n  n,j 


n+ j  -  2 


Using  the  change  of  variable  y  ■  s  -  bn,  the  integral  becomes 


“n1  4(t"  b  {.L  dv 


n  j=0 


-  o'1  T  a  n  v  ^  rn-2+j-.,n-2+3-h  ^n  h,  2  v-1  , 

"  Jo^-j  h-0  (  »  )b"  ‘sw-  bnyty*cn>  <* 


iX2  ("txxo-v  •  • 


n-2+j 


rn-2+j.h-h 


where 


xh(>')  =  /n  yh(y2  +  c  )*1  dy 


kc-O1*^*1^  log(y2  +  c_)  (h  odd) 


I  Ch - 3) / 2 ]  (-c J1  .  h  kfh  n  , 

j0  pfe*  (-«'h<5<hl)*"1^J 

l0s 

7  n 


(h  even,  cn  >  0) 


(h  even,  cr  <  0) 


[For  exanple,  I0 (y)  ■  1/*^  tan'  y/*^  if  cn  >  0; 


[1/(2/^  )]  log  ((y-/^)/(y+/^))  if  Cjl  <  0;  IjCy)  =  %  log(y  +cn); 

i2(y)  -  y  -  cn!0(y) ;  tyy)  *  Hy2  -  cnx x Cy) ;  i4(y)  -  U/3)y3  -  cni2(y);  and  so 

on.] 
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The  constants  a  ,  b  ,  and  c  are : 

n  n  n 

an  =  (n  +  l)a  +  »j(n  - 1)  (n  +  2)B  ; 
bn  =  *s{(2n  +  l)a  +  (n2-l)B)a^  ; 
cn  =  [an  -  1/4 {a  +  (n-l)g}2]a^2  . 


CNote  that  if  a  =  0,  a^  *  h(n -1)  (n  +  2)8;  b  =  (n  +  l)/(n  +  2); 
cn  =  2{(n-l)(n  +  2)6}‘1  -  (n  +  2)'2.) 

The  expressions  for  the  coefficients  d^  ^  are: 

dn  0  *  (1  +  !sn(n-l)6}[l  +  nail  -  S(x)}  +  *sn(n-l)B]  ; 

dn  .  *  -2(n2- 1) {1  +  n(n-l)B)B  -  (n  +  l)a(l  +  (3/2)n(n -1)6}{1  -  Six)} 

dn  ,  *  (n  -1)  (n  +  2)8  +  -1)  (3n2  +  6n  +  2)  [ad  -  S(  r)}  +  (n-l)B]B  ; 

"»*■ 

d  ,  *  - (n2  -1)  (n+  2)  t*s a{l  -  S(x)}  +  (n-l)B]B  ; 
d  .  =  fl/4)(n-l)2  (n  +  2)2  g2  . 
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